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Abstract 

We study the interaction between the Moon and the solar wind using a three-dimensional hybrid plasma 
solver. The proton fluxes and electromagnetical fields are presented for typical solar wind conditions with 
different magnetic field directions. We find two different wake structures for an interplanetary magnetic field 
that is perpendicular to the solar wind flow, and for one that is parallell to the flow. The wake for intermediate 
magnetic field directions will be a mix of these two extreme conditions. Several features are consistent with a 
fluid interaction, e.g., the presence of a rarefaction cone, and an increased magnetic field in the wake. There are 
however several kinetic features of the interaction. We find kinks in the magnetic field at the wake boundary. 
There are also density and magnetic field variations in the far wake, maybe from an ion beam instability related 
to the wake refill. The results are compared to observations by the WIND spacecraft during a wake crossing. The 
model magnetic field and ion velocities are in agreement with the measurements. The density and the electron 
temperature in the central wake are not as well captured by the model, probably from the lack of electron physics 
in the hybrid model. 

1 Introduction 

Bodies that lack a significant atmosphere and internal magnetic fields, such as the Moon and asteroids, can to a first 
approximation be considered passive absorbers of the solar wind. The solar wind ions and electrons directly impact 
the surface of these bodies due to the lack of atmosphere, and the interplanetary magnetic field passes through 
the obstacle relatively undisturbed because the bodies are assumed to be non-conductive. Since the solar wind is 
absorbed by the body, a wake is created behind the object. This wake is gradually filled by solar wind plasma 
downstream of the body, through thermal expansion and the resulting ambipolar electric field, along the magnetic 
field lines. For a review of the Moon-solar wind interaction, including recent findings, see Halekas et at. (2010), 
who also point out the fact that we can view the near lunar space as a plasma laboratory. Although the Moon lacks 
any significant atmosphere and ionosphere there are many interesting features to this seemingly simple solar wind 
interaction. To understand the global interaction between the solar wind and the moon we need computer models 
that can give us a three-dimensional picture that complements the in-situ plasma observations. 

Magnetohydrodynamic (MHD) modeling of the Moon-solar wind interaction has a long history, starting with 
the work of Spreiter et al. (1970), who analytically treated the special case of an interplanetary magnetic field (IMF) 
aligned with the solar wind flow direction by finding a transformation of the MHD equations into a hydrodynamic 
problem. However, there are many kinetic processes that cannot be described by fluid models, e.g., the non- 
Maxwellian particle populations in the wake region, so particle models should capture more of the physical processes. 
An approximation of the refill of the lunar wake is the general one-dimensional problem where plasma expands into 
a vacuum (Widner, 1970; Samir, 1983; Mora, 2003), and such models have specifically been applied to the refill 
of the lunar wake (Farrell, 1998; Birch and Chapman, 2001). A nice property of such models is that even if they 
are one-dimensional, they can to some degree approximate a two-dimensional model of the lunar wake, where 
time corresponds to distance downstream the wake, i.e., the one-dimensional model is applied perpendicular to the 
solar wind flow and convects with the flow. A disadvantage of one-dimensional models, e.g., along x, is that the 
magnetic field component along the x-axis is forced to be constant, from the requirement that V • B = 0. This is a 
severe restriction for modeling the lunar wake, where the field component across the wake can increase. Due to the 
computational complexity, particle in cell (PIC) models, that include electrons as particles, are for the foreseeable 
future limited to at most two spatial dimensions, so what has been investigated using such models is the interaction 
of an infinite cylinder with the solar wind (Harnett and Winglee, 2003; Kimura and Nakagawa, 2008). 
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Hybrid models, where ions are particles and electrons a fluid, are less computational demanding that full PIC 
models and are suitable for studying processes of size comparable to the ion inertial length and gyro radius. A 
two-dimensional hybrid model by Travnicek et al. (2005) focused mostly on wave activity in the lunar wake. The 
Moon-solar wind interaction is however a fully three-dimensional problem, since the IMF component perpendicular 
to the solar wind flow introduces asymmetry into an otherwise cylindrical symmetric problem. A three dimensional 
hybrid model by Kallio (2005) had a fairly coarse grid on a small region, and did not handle the wake region in 
a self consistent way (the electric field was specified in that region). Recently another three dimensional hybrid 
model has been presented in Wiehle et al. (2011). They compare observations by one of the ARTEMIS spacecrafts 
to a model run with time dependent inflow conditions, and find a good agreement. The time dependence enables 
comparison with observations even during changing solar wind conditions, but makes it more difficult to see how 
the wake forms differently depending on the solar wind conditions. 

Here we also apply a three dimensional hybrid model to the Moon-solar wind interaction, but for steady inflow 
conditions. This allow us to study how the wake depends on the direction of the IMF, and in what follows we 
describe the method and present a comprehensive global overview of the different plasma quantities in the wake. 



2 Model 

In the hybrid approximation, ions are treated as particles, and electrons as a massless fluid. In what follows we use 
SI units. The trajectory of an ion, r(t) and v(t), with charge q and mass m, is computed from the Lorentz force, 

dr dv q _ ^ 

^ =v ' -jr = -(E + vxB), 
at at m 

where E = E(r, t) is the electric field, and B = B(r, t) is the magnetic field. The electric field is given by 

E = ^- (-J 7 x B + fiQ 1 (V x B) x B - Vp e , ) (1) 

where pi is the ion charge density, J/ is the ion current, p e is the electron pressure, and po = 4tt • 10 _T [Hm _1 ] is 
the magnetic constant. We assume that the electrons are an ideal gas, then p e — n e kT ei so the pressure is directly 
related to temperature (k is Boltzmann's constant). There are several ways to handle the electron pressure (p. 
8790, Winske and Quest, 1986). Here we assume that p e is adiabatic (small collision frequency), with an adiabatic 
index, 7 = 5/3. Then the relative change in electron pressure is related to the relative change in electron density 
by p e oc |p e | 7 (Winske and Quest, 1986), and we have that 

Pe _ f n e 
PeO \n e o 

where the zero subscript denote reference values (here the solar wind values at the inflow boundary). From charge 
neutrality and p e = n e kT e we can derive that 

p e = Ap] with A = -/pT e (2) 

a constant that is evaluated using reference values of pi and T e (solar wind values). Note that 7 = 1 corresponds 
to assuming that T e is constant, and 7 = gives a constant p e . 

Finally, Faraday's law is used to advance the magnetic field in time, 

<9B ^ ^ 

^ = " VXE - 

We use a cell-centered representation of the magnetic field on a uniform grid. All spatial derivatives are dis- 
cretized using standard second order finite difference stencils. Time advancement is done by a predictor-corrector 
leapfrog method with subcycling of the field update, denoted cyclic leapfrog (CL) by Matthews (1994). An ad- 
vantage of the discretization is that the divergence of the magnetic field is zero, down to round off errors. Also, 
the discretization conserves energy very well (Holmstrom, 2011). The ion macroparticles (each representing a large 
number of real particles) are deposited on the grid by a cloud-in-cell method (linear weighting) , and interpolation of 
the fields to the particle positions are done by the corresponding linear interpolation. Initial particle positions are 
drawn from a uniform distribution, and initial particle velocities from a drifting Maxwellian distribution. Further 
details on the hybrid model can be found in Holmstrom (2010). 
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The implementation of the algorithm was done in the FLASH software framework, developed at the University 
of Chicago (Fryxell et a/., 2000), that implements a block-structured adaptive (or uniform) Cartesian grid and is 
parallelized using the Message-Passing Interface (MPI) library for communication. Further details on the FLASH 
hybrid solver can be found in Holmstrom (2011). The general solver will be included in future releases of FLASH. 

A crucial point in modeling the interaction between the Moon and the solar wind is how to choose the inner 
boundary conditions, at the surface of the moon. Here we model the Moon as an absorber of the solar wind. 
Ideally we should then remove all ions that hit the surface of the Moon, but that would introduce a discontinuity 
in the charge density, from solar wind values outside the sphere, to zero inside the sphere. However, numerical 
differentiation of discontinuous functions on a grid will lead to oscillations and instabilities in the fields. Therefore 
we choose to keep also ions inside the Moon, but we gradually reduce their weight by a factor / 5 S after each time 
step while they are inside the Moon (the mass and charge of each particle are reduced). Thus, the Moon acts as 
a gradual sink for ions, and this approach also avoids the problem of pi = in Eq. [T] However, we do not plot 
these reduced ions in what follows, only protons that have never hit the Moon are used when computing the various 
statistics. This can be seen as a form of smoothing, done in a self consistent way. Often smoothing of fields are used 
in hybrid solvers, but it is then done in an ad-hoc way, as an extra step. The downside of our approach is that we 
introduce a numerical parameter, but experience shows that the solution will not change by much when reducing 
fobs, indicating that the solution is converging. The reduction parameter cannot be chosen too small however, since 
a steep gradient in ion density will cause oscillations in the magnetic field. We have verified numerically that this 
approach keeps the IMF constant inside the Moon, and thus that no spurious currents are present there. There are 
no inner boundary conditions on the electro-magnetic fields. 

The coordinate system is centered at the Moon, with the x-axis is directed toward the Sun, so the solar wind 
flows opposite to the x-axis. The IMF is in the xy-plane, and we will from now on denote it as the IMF plane. This 
defines a plane of symmetry, since the only asymmetry in the problem stems from a non-zero component of the 
IMF perpendicular to the x-axis. The solar wind consists only of protons (and the massless electron fluid), and we 
have a rectangular simulation domain divided into cells. The -fx-face of the domain is an inflow boundary where 
we have a layer of shadow cells outside the domain that are filled with macroparticles drawn from the solar wind 
proton velocity distribution, and where the solar wind magnetic field is set. The — x-face is an outflow boundary 
where particles exiting the domain are removed, and where the magnetic field is extrapolated from the interior. The 
four other faces have periodic boundary conditions for fields and particles. The simulation domain is divided into 
a Cartesian grid with cubic cells of size 160 km. The geometry of the simulation domain is shown in Fig. [I] The 




Figure 1: The proton number density for an IMF along the direction (1, 1,0). This also illustrates the coordinate 
system and the simulation geometry. The cuts are the planes x = —10000 km, x = —30000 km, y = 0, and z = 0. 
These cutting planes are the same that later are used in Fig. 2-4. 
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extent of the simulation box is [-32.40,4.08] x [-8.32,8.32] x [-8.32,8.32] • 10 3 km w [-18,2.4] x [-4.8, 4.8] 2 i? L , 
with the Moon as a sphere of radius Rl — 1730 km at the origin. The total number of simulation macroparticles is 
about 300 million, with 121 particles per cell on average. Here we have used values of 7 = 5/3, / 6 S = 0.965, a time 
step of 0.08 s, and a final time of 170 s, when the solution has reached a steady state. The plots in what follows 
show the solution at the final time. Five subcycles of the CL algorithm are used here when updating the magnetic 
field, and the execution time is about 13 hours on 384 CPU cores. 

3 Results 

In what follows, we have used typical solar wind conditions at 1 AU (Table 4.1, Kivelson and Russel, 1995), unless 
otherwise noted. The solar wind velocity is 450 km/s, the number density is 7.1 cm -3 , and the ion temperature is 
1.2 • 10 5 K. The electron temperature is 1.4 • 10 5 K, and the IMF is in the x?/-plane with a magnitude of 7 nT. For 
these plasma parameters, the ion inertial length is 85 km, the Alfven velocity is 58 km/s, the thermal proton gyro 
radius is 66 km (the gyro time is 9 s), and the ion plasma beta is 0.6. 

3.1 Overview 

In Fig. [2] the number density (three first rows) and the ion velocities (three last rows) are shown, each for three 
different IMF directions (at 0, 45 and 90 degrees to the x-axis, in the x?/-plane). The most obvious feature of the 
interaction between the Moon and the solar wind is that a wake is formed behind the Moon, visible as a region 
of low density. We also see how the refill of the wake depends on the IMF direction. The wake refills along the 
magnetic field lines, making the wake shorter when the component of the IMF that is perpendicular to the solar 
wind is larger (Fig.|2jl). When the IMF is anti-parallel to the solar wind flow there is hardly any refill of the wake 
over the length of the simulation domain, also the flow then is cylindrical symmetric around the Moon-Sun axis. 
That the wake refill occur along the magnetic field lines is also visible in the cuts of the far wake (Fig. [2jld), with 
the refill region as a rectangle aligned with the IMF. 

Another notable feature is the density rarefaction cone that emanate from the terminator region. This region 
of low density propagate mostly above the poles in the case of a perpendicular IMF (Fig. |2]l), indicating that it 
propagates perpendicular to the magnetic field lines. This will be discussed in more detail when we look at the 
associated reduction in magnetic field strength (Fig. |3|. 

We now examine the ion velocities, computed from the ion current and the charge density as J//p/, as is usual 
for particle in cell models, since the ion current density and the ion charge density are quantities that are available 
on the computational grid. The small white region downstream of the Moon is a region without any particles that 
has not hit the Moon, thus no velocities can be computed there. The x-velocity is shown as the deviation from 
the solar wind velocity, when the IMF is perpendicular to the flow, in Fig. |2j4. We see a slight slow down of the 
solar wind in the rarefaction cone. There is also a slow down behind the +z hemisphere (and a speed up behind 
the — z hemisphere) seen in Fig. [5]4b. This is probably a geometric effect, where part of the velocity distribution 
function is emptied by collisions with the lunar surface. The y- velocity (Fig.[2j5) shows how the wake symmetrically 
refills from both sides of the wake. The ions move inward not only in the central part of the wake, but also in the 
rarefaction cone. For the case of an aligned IMF, the z- velocity (Fig. [2]6) also shows a flow in the rarefaction cone 
toward the wake center. There we also see a flow of ions at the wake boundary in the IMF plane, and there seem 
to be a small kinetic effect of alternating upward and downward going ions in the upper and lower part of the wake 
(Fig.[2]6b). This seems associated with the wave pattern in density inside the rarefaction cone (Fig. [2] 1-3), possibly 
caused by oscillations at the wake boundary. 

The magnetic field magnitude (three first rows) is shown for three different IMF directions in Fig. |3j where also 
each of the vector components (rows 4-9) are shown for two different IMF directions (perpendicular and aligned to 
the solar wind flow). 

The most prominent feature is the enhanced magnetic field in the wake. For a perpendicular IMF (Fig.|3]l) this 
region spreads out in the IMF plane, in the same region where we in the density saw the wake refill. For the aligned 
IMF the magnetic field increase is larger, and confined to the narrow wake (Fig. [3J3). There are also structured 
fluctuations in the far wake region in the IMF plane. Then there is a cone of lower field strength in the case of an 
aligned IMF (Fig.[3j3), corresponding to the cone of reduced density (Fig.|2]3). For the perpendicular IMF the cone 
is partial, and only extends perpendicular to the IMF, toward +z and —z (Fig. [5Jl). Between the central region 
of enhanced field and the cone of lower field, there are field fluctuations parallel to the cone, as is also seen in the 
density (Fig. [2] lb, 2b, 3b). Examining these density plots, the waves seem to originate at the wake boundary and 
could be caused by oscillations at this interface between high and low ion density. 
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Figure 2: Proton number density and velocities in different planes, for different upstream IMF conditions. Arrows 
show the direction of the IMF. Top three rows show number density for different upstream IMF directions. Along 
the directions (x,y) = (0, 1), (1, 1), and (1,0), respectively. The geometry is repeated in the lower three rows (4-6) 
for the proton velocities u x — u swi u y , and u z . The columns are from the left, cuts in the planes z = (seen from +z), 
y = (seen from x = —10000 km, and x = —30000 km (seen from -\-x). Fig. 1 provides a three-dimensional 
illustration of the cutting planes. The two vertical lines in l-3b show the position of the cuts perpendicular to the 
x-axis. 



The magnetic field z-component for the case of a perpendicular IMF (Fig. [3] 8) exhibits a pinching of the field, 
consistent with earlier predictions (Fig. 3, Owen, 1996) and model observations (Kallio, 2005). The field is bent 
toward the IMF plane in the wake (Fig. [3j8c,8d). However, the x-component of the magnetic field is also perturbed 
(Fig. [3j4). Close to the IMF plane, the field is bent toward the moon on one side of the wake, and away from the 
moon on the other side of the wake (Fig.[3j4a). Note also that these bends in the magnetic field are not confined to 
central parts of the wake, but expand along with the rarefaction in magnetic field (Fig.[3jld) and density (Fig. [2] Id). 
Thus, for a perpendicular IMF the magnetic field in the wake is not only bent toward the IMF plane, but also toward 
the moon, as shown in Fig. |1J 

The Mach cone of decreased density and magnetic field makes an angle to the solar wind flow direction of about 
10 degrees, as seen in Fig. [2] and Fig. [3J Since the outer edge of the cone is diffuse, it is impossible to give a 
more exact angle. The ion-acoustic wave velocity is v$ ~ ^2^kBT e /rrii « 58 km/s, and the Alfven velocity is 
also va ~ 58 km/s for these solar wind conditions, corresponding to angles of 7 degrees. We see that the density 
depletion travel perpendicular to the magnetic field. This together with the fact that the propagation velocity is 
larger than both the ion- acoustic and the Alfven velocity suggests that the fast magnetosonic wave, with a phase 
velocity of \Jv\ + v 2 s leading to an angle of 10 degrees, is responsible for the creation of the Mach cone in the hybrid 
model. This is consistent with the observation of Wiehle et al. (2011) that all three MHD modes (fast, Alfven and 
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slow) will be manifested in the lunar wake. 

In Fig. [5] the electric field magnitude (top row) is shown along with the three vector components (three lower 
rows) for the case of an IMF perpendicular to the solar wind flow direction. What is shown is the electric field 
with the solar wind convective electric field, E sw = — u sw x B sw , (directed along the z-axis) subtracted. We do 
not show the electric field for the case of an IMF aligned with the x-axis, since it is indistinguishable from zero 
using the same color scale. We can note that the regions of enhanced electric field magnitude coincide with those of 
enhanced magnetic field magnitude shown in Fig. [3}l. Also, there are similarities in the components of the electric 
and magnetic fields. E y is very similar to B z (Fig7[3]8), and correspondingly E x is similar to B y (Fig. |3]8). 

3.2 Comparison with WIND observations 

We here compare the hybrid model results with observations by the WIND spacecraft on December 27, 1994, when 
it traversed the lunar wake approximately perpendicular to the solar wind flow, in the plane of the IMF, at a 
distance of x = —6.5 Rl ~ —11300 km. Ogilvie et at. (Fig. 1, 1996) plots several plasma parameters during the 
wake crossing. Apparently the solar wind conditions changed during the observation. This is most apparent in 
the plasma number density that decreased by more than a factor of two from the undisturbed solar wind on the 
inbound trajectory to the solar wind on the outbound trajectory. Therefore we make two runs of the hybrid model 
for the two sets of solar wind conditions, here denoted before and after. The solar wind conditions used for the 
before (after) case are as follows. A solar wind velocity of 470 (500) km/s, a number density of 5 (2) cm -3 , an ion 
temperature of 0.9 (0.5) • 10 5 K, and an electron temperature of 1.5 (2) • 10 5 K. The IMF was (6.39, 2.33, 0) nT for 
the before case, and (5.57,5.02,0) nT for the after case. The time step was 0.08 (0.05) s and f b s = 0.98(0.985). 
Since the plasma parameters are different in the two cases, these numerical parameters have different values to 
ensure the stability of the solutions. 

In Fig. [6] we show the model magnetic field magnitudes, compared with the observation by WIND. The observed 
magnetic field magnitude matches well to the model values, if we assume that there was a change in the upstream 
solar wind conditions around wake entry (there we switch from comparing the before model results to comparing 
the after model results with the observation). Before the wake crossing there is a drop in the observed magnetic 
field strength that is also seen in the model. This is the rarefaction wave that we saw in Fig. [3] earlier. The IMF 
in the before case has the largest component along the x-axis and therefore the magnetic field magnitude should 
best match Fig. [3j3a, with a pronounced rarefaction cone. Then there is an increase in the central wake which is 
reproduced in the after model. Followed by a gradual decrease at wake exit, but no rarefaction wave. For the after 
case the IMF is at about 45° away from the x-axis and thus correspond to Fig.[3j2a, with a much weaker rarefaction 
cone. In the observation there are some large variations further away from the wake that are not present in the 
model. This could be upstream solar wind disturbances. 

In Fig. [7] the proton velocity observed by WIND is compared to the velocity in the two hybrid model runs. 
Again we have fairly good agreement between the model and observations, assuming a shift in solar wind conditions 
around wake entry. There are more variations in the observed velocities, but the locations of the velocity changes, 
and the slopes match the model. In the central wake there is however no data from the model since the model 
density is too low. This can be seen when we compare the proton number density from the model with observations 
in Fig. [8| We see that the general features of the observed density profile across the wake is found in the model, 
including a shift of the minimum density toward the outbound part of the trajectory, if we assume the shift in 
solar wind conditions around wake entry. However, the drop off in density is quicker in the hybrid model, and is 
much below the observations in the central part of the wake. The cause of this could be the incomplete treatment 
of electron dynamics in a hybrid model. Since the electrons are a massless charge neutralizing fluid, there are no 
forces from the electrons on the ions, except for the electron pressure gradient term in the definition of the electric 
field (Eq. [I]). In reality, the refill of the wake is governed by the electrons. Due to their higher thermal velocity they 
will speed ahead of the ions to refill the wake, setting up an ambipolar electric field due to the charge separation 
that will accelerate the ions into the wake (Halekas et al., 2010). 

Similarly, the lack of electron dynamics in the hybrid model will affect the electron temperature, T e , in the wake. 
Observations show an increased electron temperature in the wake (Fig. 1, Ogilvie et al, 1996). But from Eq. [2J 

1 2 /3 

and the ideal gas law, T e ~ p]~ = p/ for 7 = 5/3, so T e in the hybrid model is then lower in the wake. 

4 Discussion 

Since the only deviation from axial symmetry is the component of the IMF that is perpendicular to the solar wind 
flow, we can see that the two extreme cases are when the IMF is perpendicular to the flow, and when the IMF is 
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aligned with the flow. For a constant IMF magnitude, the flow and field for cases when the IMF is at an angle to 
the flow will be a mixture of these extreme cases, as seen in Fig. [2j2 and Fig. [3j2. 

There are many processes related to electron dynamics that are not captured by a hybrid model. The process 
of plasma expansion into a vacuum has a lot of details, and instabilities that can form (Birch and Chapman, 2001). 
Also, the region near the lunar surface probably has strong electric fields, and charge separation effects (Farrell et 
al, 2008). However using the hybrid approximation it is possible to model the global three-dimensional interaction 
between the Moon and the solar wind, so full PIC models and hybrid models are complementary tools in the study 
of this interaction. Also, the Moon-solar wind interaction, although seemingly simple, is complex enough to be a 
good test of the applicability of different models; fluid, hybrid, and PIC models. 

We have here considered the simplest possible model of the Moon in the solar wind — a spherical sink of solar 
wind protons. There are many physical processes that make the real interaction more complex than that. Two 
examples are magnetic anomalies that deflect the solar wind and will cause downstream disturbances (Lue et al., 
2011), and the reflection of solar wind protons by the dayside lunar surface (Holmstrom et al, 2010). 

5 Summary and conclusions 

We have presented a fully self consistent three-dimensional hybrid plasma model of the interaction between the 
Moon and the solar wind and shown the resulting global ion fluxes and fields for typical solar wind conditions and 
different IMF directions. Magnetic field enhancement in the wake, and a rarefaction wave in density and magnetic 
field are seen in the model. The two extreme cases are that of an IMF perpendicular to the solar wind flow, and 
that of an IMF that is parallel to the solar wind flow. We find that the wake for intermediate IMF directions can 
be seen as a superposition of these two states. Some kinetic effects were also observed in the model. There are 
density and field variation in the far wake, maybe caused by an ion beam instability from the refill of the wake. 
Visible are also kinks in the magnetic field at the inner boundary of the rarefaction cone for an IMF perpendicular 
to the solar wind flow direction. 

Observations by the WIND spacecraft were fairly well reproduced by the hybrid model, showing that a hybrid 
model with the moon as a spherical sink of ions is able to represent many global features of the interaction. However, 
the model central wake density and electron temperature were different from observations, probably due to the lack 
of electron dynamics in the hybrid model. 
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Figure 3: The magnetic field magnitude, and components in different planes, for different upstream IMF conditions. 
Arrows show the direction of the IMF. Top three rows show magnetic field magnitude for different upstream IMF 
directions. Then the magnetic field x-component for an IMF perpendicular to and opposite the solar wind flow is 
shown on row 4 and 5. Similarly, the ^-component is shown on row 6 and 7, and the z-component on row 8 and 9. 
The geometry of the cuts in the different columns are the same as in Fig. 2. The colorbar is for the row immediately 
to the left, and for all following rows, until the next colorbar. 
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Figure 4: Magnetic field lines 10000 km behind the Moon for an IMF along the y-ax.is (same as in Fig. [3| 1) . (a) 
seen along the +a>axis, and (b) along the -z-axis. The color scale shows the magnitude of the magnetic field 
similarly to Fig. [3] To enhance the bending of the field the magnetic field ^-component was multiplied by 4, and 
the ^-component by 2, before drawing the field lines. 




Figure 5: The electric field magnitude, and the three components in different planes, for an IMF that is perpendicular 
to the solar wind flow direction. Arrows show the direction of the IMF. For the magnitude and ^-component, the 
solar wind convective electric field, E sw = — u sw x B sw , has been subtracted. The cuts in the different columns are 
the same as in Fig. 2 and 3. 
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Figure 6: The hybrid model magnetic field magnitude [nT] along y at x = —6.5 Rl =-11245 km for two different 
upstream solar wind conditions (labeled before and after), compared to WIND observations from Ogilvie et al. (Fig. 
1, 1996). The vertical lines show the location of the optical shadow, and the x-axis show time on December 27, 
1994, as in the original plot. 




Figure 7: The hybrid model proton x- velocity magnitude [km/s] for two different upstream solar wind conditions 
(labeled before and after), compared to WIND observations. Otherwise similar to Fig. [6] The reason that the 
WIND velocity is double valued in the center of the wake is that the velocity of the two proton beams refilling the 
wake from each side is computed separately, as explained in the original publication (Ogilvie et al., 1996) 
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Figure 8: The hybrid model proton number density [cm -3 ] for two different upstream solar wind conditions, 
compared to WIND observations. Otherwise similar to Fig. [6] 
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